#include "../include/Element.hpp"
#include "../include/Region.hpp"
#include "../include/Grid.hpp"
#include "../../input/parameters.hpp"
#include "../utils/load_point.hpp"
#include "../utils/load_region.hpp"
#include "../utils/load_element.hpp"

int main()
{
    Region ref = load_region();

    Point p0(0, {1.0, 1.0}, 1);
    Point p1(1, {2.0, 1.0}, 1);
    Point p2(2, {1.0, 2.0}, 1);

    Element e0(0, {p0, p1, p2}, 1);

    std::function<double(double, double)> f = [](double x, double y)
    {
        return 1;
    };
    double rst = e0.integral(f, ref);
    std::cerr << rst << std::endl;
    return 0;
}
